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ABSTRACT 



Following renewed interest, the problem of whether the cosmological expansion af- 
fects the dynamics of local systems is reconsidered. The cosmological correction to the 
equations of motion in the locally inertial Fermi normal frame (the relevant frame for 
astronomical observations) is computed. The evolution equations for the cosmological 
perturbation of the two-body problem are solved in this frame. The effect on the orbit 
is insignificant as are the effects on the galactic and galactic-cluster scales. To appear 



in the Astrophysical Journal 
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1 Introduction 



A recurring issue in cosmology concerns the nature and extent of the cosmological expan- 
sion. If expansion were to occur in proportion and in every minutia of detail to every 
element of the universe, then every clock and every measuring apparatus of distance 
would be altered in proportion. If in addition, the laws of physics were to remain unal- 
tered in the process, the very concept of expansion would lose its meaning as it would be 
intrinsically unobservable. This is not what is being contemplated for our universe since 
we observe a systematic redshift of distant galaxies and hence we are able to deduce that 
there is an expansion in progress, at least on the cosmological scale. While the effect 
is actually registered in the small distances of a wavelength of light, this is simply an 
imprint of the expansion at the largest (Hubble) scale. 

Recently (Anderson 1995; Bonnor 1996) there has been a revival of interest in the 
question as to whether the cosmological expansion also proceeds at smaller scales. There 
is a tendency to reject such an extrapolation by confusing it with the intrinsically un- 
observable "expansion" (let us refer to this as "pseudo-expansion") described above. 
By contrast, the metric of Friedman-Robertson- Walker (FRW) in general relativity is 
intrinsically dynamic with the increase (decrease) of proper distances correlated with 
red-shift (blue-shift). It does so on any scale provided the light travel time is much 
longer than the wave period. Thus, the cosmological metric alone does not dictate a 
scale for expansion and in principle, it could be present at the smallest practical scale 
as real- as opposed to pseudo-expansion, and observable in principle. 

However, it is reasonable to pose the question as to whether there is a cut-off at which 
systems below this scale do not partake of the expansion. It would appear that one would 
be hard put to justify a particular scale for the onset of expansion. Thus, in this debate, 
we are in agreement with Anderson (1995) that it is most reasonable to assume that the 
expansion does indeed proceed at all scales. However, there is a certain ironical quality 
attached to the debate in the sense that even if the expansion does actually occur at all 
scales, we will show that the effects of the cosmological expansion on smaller spatial and 
temporal scales would be undetectable in general in the foreseeable future and hence 
one could just as comfortably hold the view that the expansion occurs strictly on the 
cosmological scale. 

The question of whether the expansion of the universe affects local systems like 
clusters of galaxies or planetary systems was first raised many years ago and has received 
continued scrutiny (McVittie 1933; Jarnefelt 1940, 1942; Pachner 1963; Dicke & Peebles 
1964; Callan et al. 1965; Irvine 1965; Noerdlinger & Petrosian 1971) with the most 
recent consideration by Anderson (1995) who extends the question to the stellar scale 
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and even below this. The recurrent attention paid to this issue indicates that to this 
point a definitive answer is still lacking. However, it is our sense that the prevalent 
perception is that the physics of systems which are small compared to the radius of 
curvature of the cosmological background is essentially unaffected by the expansion of 
the universe. 

In the presence of spherical symmetry, the analysis of a spherical cavity embedded 
in an FRW universe is well known: as a consequence of Birkhoff 's theorem, the metric 
inside the spherical cavity is the Minkowski one, and the physics is the same as in flat 
space (Einstein & Straus 1945; Schiicking 1954; Dicke & Peebles 1964; Callan et al. 1965; 
Bonnor 1996). However, when spherical symmetry is absent, a satisfactory quantitative 
answer is missing in the literature, and certain statements about small systems being 
sheltered from the cosmological expansion recur (see e.g. Misner et al. 1973 and the 
discussion in Anderson 1995). Noerdlinger & Petrosian provide a quantitative treatment, 
but it is limited to the particular problem of the collapse of galaxy clusters. Anderson's 
(1995) paper employs the Einstein-Infeld-Hoffmann method to derive the cosmological 
corrections to the equations of motion of a system of particles subject to external forces. 
When the dynamics of a single particle are considered, the correction to the particle's 
acceleration is found to be proportional to the velocity of the particle. However, it is not 
clear how to relate the coordinates used by Anderson to the coordinates used by observers 
making astronomical observations. This is an important issue because the computation 
does not provide a coordinate-invariant quantity, but rather a correction to the 3- 
dimensional equations of motion, that are dependent upon the chosen coordinate system. 
Bonnor (1996) studies a distribution of pressureless charged dust in equilibrium between 
electrical repulsion and gravitational attraction, and concludes that it participates in the 
universal expansion. 

A qualitative answer to the problem of whether local systems are affected by the 
expansion of the universe is easily provided if one considers the equivalence principle 
and its geometric formulation. Although the cosmological expansion is described by the 
time-dependent scale factor in the FRW metric, and we believe affects lengths at all 
scales, the curved spacetime manifold can be locally approximated by its (flat) tangent 
space at every spacetime point p. This approximation is valid only in a neighborhood 
U(p) of the point p considered; the error involved in the approximation increases with 
the size of the neighborhood U {p) , and the approximation breaks down completely when 
the size of U (p) becomes comparable to the radius of curvature of spacetime (the Hubble 
radius in the case of an FRW spacetime). From the physical point of view, the tangent 
space at p describes the spacetime seen by a freely falling observer in the so-called locally 
inertial frame (hereafter called "LIF"- see Landau & Lifshitz 1989). This frame is the one 
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in which astronomical observations are carried out. Thus, the effect of the cosmo logical 
expansion is seen to be negligible locally and grows in significance with distance, reaching 
full import on the cosmological scale. This conclusion is qualitative, and is certainly 
well-known to most relativists but, to the best of the authors' knowledge, has yet to be 
well-formulated quantitatively. In earlier treatments, the coordinate systems adopted 
do not correspond to those used by a physical observer. 

The purpose of the present paper is to provide a clear quantitative answer to the 
problem. The motion of a particle subject to external forces in the (approximate) LIF 
using Fermi normal coordinates is analyzed. It is the locally inertial frame based on 
a geodesic observer and it continues to be locally inertial following the observer in 
time. This is the frame in which astronomical observations are performed, and we 
compute the corrections to the dynamics due to cosmology. In this paper, we assume that 
homogeneous isotropic expansion is actually universal and we analyze the consequences 
of this assumption. 

The plan of the paper is as follows: in Sec. 2 the equation of timelike geodesies in the 
LIF in an Einstein-de Sitter universe is investigated and the cosmological perturbations 
to the 3-dimensional equations of motion in the LIF are derived. In Sec. 3 the orders 
of magnitude of the effects in realistic astrophysical systems are estimated, and it is 
demonstrated that they are very small and unobservable with present and foreseeable 
technology. In Sec. 4 the two-body problem in the LIF is studied in detail using the 
correction to the equations of motion computed in Sec. 3, thus providing a solution to 
the evolution equations for the perturbations. It is shown that cumulative effects of the 
cosmological expansion on the present orbital radius of the earth and its orbital motion 
are essentially negligible. Section 5 contains a discussion and the conclusions. 



2 Equations of motion in the LIF 



In this section we find the equations of motion for a particle in the LIF using the geodesic 
deviation equation. We refer the reader to the Appendix for the details of the calculation. 
The metric in FRW coordinates for an Einstein - de Sitter universe is given by 



ds 1 



-dt 2 + a{t) 2 dx 2 + dy 2 + dz z 



(2.1) 



where a(t) is the scale factor. Consider an observer whose world line is the geodesic 
t — t and r — 0. In Fermi normal coordinates, the metric on the geodesic is g^ u = r)^ u 
and a parallely propagated orthonormal tetrad is given by 



;i, 0,0,0) 



(0,1,0,0) 



(2.2) 
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ef = (0,0,1,0) ef = (0,0,0,1). (2.3) 
The FRW basis vectors e^ RW are related to the via 



< = A/e™ ( 2 - 4 ) 

where =diagonal(l,a _1 , a -1 , a -1 ). The Riemann tensor in Fermi normal coordinates 
{xp}, along the geodesic, is given by (Greek indices range from to 3, Latin indices 
from 1 to 3) 

T)F A O" A A A Ul A K TjFRW (r) r\ 

n-aPfjiv — lv a IV /3 ■'V " -"-o-Awk 

(where the superscript F denotes quantities in Fermi normal coordinates) and the 
geodesic deviation equation 

becomes 

+ i^ 4 = (2.7) 

in Fermi normal coordinates, since T^J^ = 0. Thus, to lowest order in x M and dx^/dr 
(the order to which the geodesic deviation equation is valid) , the equations of motion in 
Fermi normal coordinates are (see the Appendix) 

d x j? ( ci\ 



where a overdot denotes differentiation with respect to the comoving time t. 



3 Order of magnitude estimates 

In this section, the order of magnitude of the effect created by the cosmic expansion 
on the dynamics of local systems is estimated. Astronomical systems for which the 
velocities involved are non-relativistic are considered. The present value of the age of 
the universe is taken to be 6.3 x 10 17 seconds. 

Acceleration on the scale of the Solar System 
It is sufficient for our purposes to use the present value of the average size of the earth-sun 
system, i.e. the astronomical unit tq of 1.5 x 10 11 m and the present orbital frequency 
ujq of 2 x 10 -7 sec -1 . From eq. (|4.10|) of the following section, the correction to the 
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acceleration for this distance and frequency due to the cosmological expansion at the 
present matter-dominated (a(t) oc t 2 / 3 ) epoch is 

5f = = -3.17 x 10- 47 m/sec 2 . (3.1) 

This is to be compared to the predominant gravitational acceleration of the earth towards 
the sun 

g = = 6 x l(T 3 m/sec 2 (3.2) 

r o 

which completely overwhelms the effect of the cosmological expansion by 44 orders of 
magnitude. This is in qualitative agreement with a cruder order of magnitude estimate 
in Lightman et al. (1975). 

Acceleration on the galactic scale 
As an example, consider our galaxy, a spiral in which the sun is located at r = 8.5±1 Kpc 
from the center and has orbital velocity v = 220±15 Km s" 1 (Binney & Tremaine 1987). 

Thus, the gravitational acceleration of the sun towards the center of the galaxy 
g = Vq/tq ~ 1.9 x 10~ 10 m/sec 2 . Since the orbital period is ~ 2 x 10 8 yrs, the angular 
velocity u is of the order 1CT 15 sec -1 . Thus, the correction of the acceleration of the 
sun due to the cosmological expansion 

r = -3^1 (3 - 3) 

is of the order 10~ 21 m/sec 2 at the present epoch which is 11 orders of magnitude smaller 
than the galactic g. 

Acceleration on the galactic cluster scale 

Assuming the core radius of a galaxy cluster (r ~ 250 Kpc) and the line of sight velocity 
dispersion a ~ 800 Km s" 1 (Binney & Tremaine 1987) and assuming that the galaxies at 
the edge of the core of the cluster are in orbit around the centre of the core with velocity 
a, they are subject to the gravitational acceleration g = (fo) 2 /To — 8 x 10~ n m s~ 2 . 
This is to be compared with the correction due to the cosmological expansion 

f = ---^2 = -5.6 x 10" 18 m/sec 2 . (3.4) 

Ob 1) 

While this is 7 orders of magnitude smaller than the galactic cluster g and thus of 
considerably greater relative significance than was found for the galactic and the solar 
system scales, it is still nevertheless essentially ignorable. 
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4 Cosmological corrections to the two— body prob- 
lem in the LIF 



The effects of the expansion of the universe on the dynamics of local systems are exem- 
plified by the corrections induced in the two-body problem. The two-body problem in 
a cosmological background has been analyzed in previous papers (McVittie 1933; Dicke 
& Peebles 1964; Noerdlinger & Petrosian 1971; Anderson 1995) with differing results. 
McVittie (1933) reached the conclusion that the orbital radius stays constant for an 
observer using coordinates fixed in the solar systemF]. Dicke & Peebles (1964) used a 
conformal technique to show that the coordinate radius of the orbit decreases as the 
inverse of the scale factor and that the proper radius stays constant. Noerdlinger & 
Petrosian (1971) considered the two-body problem inside a cluster and found that, in a 
dust -dominated FRW universe, the time derivative of the average orbital radius obeys 
the equation 

(r) = H(r) , (4.1) 

W 1 + 4e W ' K J 

where e is the ratio of average energy densities in the cluster and in the rest of the 

universe (compare eq. (8) of Noerdlinger & Petrosian 1971). According to this result, the 

orbital radius increases and the effect is proportional to \Hqx 1 \. The more recent result of 

Anderson (1995) agrees with that of Dicke & Peebles (1964), but in different coordinates. 

A comparison of all these results is rendered difficult by the different coordinate systems 

adopted in the different studies. Moreover, no treatment of the problem was given in the 

LIF, which is the frame of reference relevant for astronomical observations performed by 

a freely falling observer. In fact, the 3-dimensional equations of motion of a particle are 

not coordinate-invariant and, like the equations of motion themselves, the correction 

due to the cosmic expansion is dependent upon the frame employed. In this section, 

we apply the results obtained in Sec. 2 to compute the perturbations of the two-body 

problem in the LIF in an expanding, matter-dominated Einstein-de Sitter universe. For 

simplicity, we restrict ourselves to the case of circular orbits, in which the equation of 

motion for the two-body problem takes the form 

d 2 r a _ GM , . 

^-- f= 4 - 2 

where M is the mass of the central object. In this section we only use quantities defined 
in the LIF, and we drop the subscripts. Cylindrical coordinates (r, 9, z) are used, with 

1 See Ferraris et al. (1996) for a modern criticism of McVittie's coordinates in astrophysical 
applications. 
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associated unit vectors e^, eg and e^. Since the perturbation of the central force is also 
central, the motion is again confined to the unperturbed orbital plane. We consider the 
perturbation of the orbital coordinates r, 9 given by 



r(t) = r + Sr(t) , 
9(t) = oo t + 59(t) . 



Substitution into eq. ( |4.2|) yields 



r 59 + 28ruiQ 
2u 







59 



5r 



r 



5r 



3ujQSr — 2coor S9 r = 



where 



a 
a 



2 

9t2 



in a matter-dominated universe. Combining eqs. ( |4.6| )- 



yields 



5r 



u Q 8r + — r 



0. 



(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 

(4.8) 
(4.9) 



It is easy to show that Sr is negligible relative to the other two terms in eq. ( |4.9| ) when 
t is of the order of the age of the universe and luq 1 is of the order of a year. Thus, we 
find that 

2 

1 



r{t) ~ To 



9t 2 



(4.10) 



For the earth-sun system, we take tq ~ 1.5 x 10 11 m (although strictly speaking, it is 
actually the value of r as t approaches infinity). The value of t is taken as the age of 
the universe, approximately 2 x 10 10 years or 6.3 x 10 17 sec. The angular frequency u 
is taken to be for the earth-year, approximately 2 x 10~ 7 sec -1 . 
From eqs. (|4~6| ) and (|4.10|) , 

59=^—. (4.11) 

Thus: 

1) The angular velocity decreases with time: as t approaches infinity, to approaches 
ujq so strictly speaking, u is actually the terminal angular velocity. 
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2) The orbit size grows with time: as t approaches infinity, r approaches ro- 
Now consider the fractional rate of change of frequency: 



= -8.9 x lO^sec -1 = -2.8 x lO^yr" 1 • (4.12) 

present 



This can be compared to the observed rate of variation of the orbital period of the moon 
about the earth, 2.22 ±0.35 x 10~ n yr _1 , which is larger by approximately 22 orders of 
magnitude. The cosmological effect is not significantly different at the birth of the solar 
system. For t, we would use the present time minus the age of the solar system which 
is still of order 10 17 sec. Thus the rate at birth of the solar system was not significantly 
different. 

For the fractional change in radius of the orbit, we use eq. (4.6) to find 

Z--f, (4.13) 

which gives the same kind of insignificant rate of radius growth with the expansion of 
the universe. Over the life span of the solar system, of order 10 17 sec, the fractional 
change in radius was a mere 10 -24 . 



5 Discussion and conclusions 

The effect of the cosmic expansion on the dynamics of local spherically symmetric sys- 
tems is well-known (Einstein & Straus 1945; Dicke & Peebles 1964; Callan et al. 1965; 
Bonnor 1996). In the non-spherical case, it is generally recognized that the expansion 
of the universe does not have observable effects on local physics, but few discussions 
of this problem in the literature have gone beyond qualitative statements. A serious 
problem is that these studies were carried out in coordinate systems that are not easily 
comparable with the frames used for astronomical observations and thus obscure the 
physical meaning of the computations. Moreover, different treatments lead to appar- 
ently conflicting results, as in the case of the two-body problem. This is the reason why 
the computations of Sees. 2 and 4, performed in the LIF, are particularly relevant to the 
problem. While it is reasonable to assume that the time dependence of the scale factor 
in the FRW metric ( j2 . 1| ) affects lengths at all scales in principle (see the discussion in 
Anderson 1995; Bonnor 1996), the magnitude of the effect in the LIF is the physically 
relevant one, and its computation constitutes the essential aspect of this work. 
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The computation of the cosmological correction to the local equations of motion 
performed in Sec. 2 allows one to estimate numerically the magnitude of the correction 
to the acceleration of a particle subject to external forces. The numerical estimates 
obtained in Sec. 3 suggest that the correction is extremely small and unobservable for 
galaxy clusters, galaxies and the solar system, and negligible for smaller systems such 
as stars and even more so for molecules and atoms (cf. Anderson 1995). When the 
cosmological correction to the local equations of motion is applied to the Newtonian 
two-body problem, the evolution equations for the perturbation of the orbit can be 
solved. It is found that the cumulative effect of cosmological expansion on the radius and 
angular motion of the sun-earth system is also negligible. The cosmic expansion plays 
an increasingly important role for systems whose sizes and lifetimes become increasingly 
comparable to the Hubble radius and to Hubble times respectively. In this case, the 
approximation used in this paper becomes invalid. It is well-known that the cosmological 
expansion must be taken into account, for example, in the fluid dynamical treatment 
of the formation of structures in the universe (Weinberg 1972). As a conclusion, it 
is reasonable to assume that the expansion of the universe affects all scales, but the 
magnitude of the effect is essentially negligible for local systems, even at the scale of 
galactic clusters. 
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Appendix 



In this appendix we find the transformations from FRW coordinates to Fermi normal 
coordinates. We also find the metric to order \x\ 2 and the equations of motion to lowest 
order. 

Consider an observer whose world line is the geodesic r = 0. To find the Fermi 
normal coordinates of a point P = (tFRW-> x frw) we find the unique spacelike geodesic 
which goes through the point P and intersects r = orthogonally. For a sufficiently 
small region about r = 0, such a unique geodesic is guaranteed to exist. Let Q be the 
point of intersection between r = and this geodesic and let the geodesic parameter r 
be zero at Q. The initial velocity vector T M = dx^/dr\ T=0 is chosen so that the geodesic 
reaches P at r = 1. The Fermi normal time tp is taken to be the proper time from 
the initial cosmo logical singularity to the point Q along the geodesic r = 0. The Fermi 
normal spatial coordinates are given by the projection of T M onto the orthonormal triad 
ef a) , where efoefc, and ef 3) 

The geodesic equation has solutions of the form 



ef^, and point in the x, y, and z directions respectively. 



x 



x (t, c m ) 



t = t(r, c m ) 



(A.I) 



where {c m } is a set of eight constants. From the above discussion we have the following 
conditions 

x k (0,c m ) = t F = t(0,c m ), (A.2) 



dt . 



dec ^ 

xp = a(t F ) — (0,c m ) . 



(A.3) 



This set of equations allows us to solve for c m = c m (x F ,t F ). Substituting this into 

X FRW = ^(ij °m) tFRW = C m ) (A.4) 

gives the required transformations 



X FRW ~ X FRw{ x1 Fi^f) 



t 



FRW — tFRW 



{X™, tp) 



(A.5) 



The geodesic equations are 
dH 



dr 2 



+ aa 



dr I \dr I \dr j 



= 0, 



(A.6) 
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and 



From eq. ( |A.7|) we have 



d 2 x k a dt dx k 

+ 2- — — =0 



dr 2 



a dr dr 

dx k C\ 
dr a 2 ' 



where the C\ are constants. From equation (|A.6|) we have 



dt _ 
dr~ = \ 



C- 2 



ldl s 



,2 ' 



(A.7) 
(A.8) 

(A.9) 



where C 2 is a constant. 

We now specialize to FRW spacetimes with 

a(t) = (at) n 



(A.10) 



where a and n are constants. The above differential equations have the power series 
solutions 



t(r) =t + 



\ 2 (at Q ) 



n\d\' 



[2n+l)n\C 1 \ : 



2n 



2(at ) 2 % T + 6(at ) 2 H 2 \ 2 + {at,) 2 " 



C2 + -¥*Ll* + 0(7*), 



(A.11) 



(at ) 2 " T (crfo) 2 ^^ K>) 2 " 



C 2 + 



T 



(2n + 1)C 2 + (3n + 1) 



|C*i| 2 



(crfo) 2 



3(at ) 2n *§ 



r 3 + 0(r 4 ) 



(A.12) 



where £ an d Xq are constants. Now tp = t(r = 0) gives t F = t , x k {r = 0) = 
gives Xq — 0, Xp = a(to)dx k /dr\ T= o gives C k = (at F ) n x F , and dt/dr\ T=0 = gives 
C*2 = — \x F \ 2 . Thus, using tFRw = ^(t" = 1) and x k FRW = x k {r = 1), we have 

U +o(|*jH 4 ) 



t 



FRW 



n\x F \ 
~ 1 2t f 



and 



x 



FRW 



Xl 

(atj 



1 , n 2 \x F \ 



3t 2 F 



+ 0(\x F \' 



(A.13) 
(A.14) 



12 



Note that to lowest order x k F = a(t frw ) %frw an d tp = tFRw, so that to lowest order 
Fermi normal coordinates are just "physical" coordinates in FRW spacetime. 

The spatial components of the geodesic equation, to lowest order in x k and x k , are 



which is identical to 



d 2 xp a _ 

~1T2 Xf = 

dt z a 

. The metric in Fermi normal coordinates is 



(A.15) 



ds 2 = - 



1 



n(n — 1)\xf\ 



t 2 

<F 



dtp + 



Sh\1 



n 2 \xp\ 
~3vT 



+ 



IFF 

n xi x 



3t F 



dx%dx l F . (A. 16) 



This can be written, to lowest order in Xp, as 

'4 



ds — (1 + R mm x F Xp)dt F 



^oijm x f x7 f ) dtpdxp 4 



g lh iljm J/ F u 'F 



since, to lowest order, the nonzero components of the Riemann tensor are 



r>F _ r>F 



a 
a 



n(n — 1) 

l F 



and 



R 



xyxy 



R 



yzyz 



aj 



n 

l F 



d% p< doc jp 
(A.17) 

(A.18) 
(A.19) 



(plus components related to these by symmetry). This expression is identical to the 
metric in Fermi normal coordinates given by Manasse and Misner (1963). 
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